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A charge moving freely in orbit around the Earth radiates according to Larmors formula. If the 
path is closed, it would constitute a perpetuum mobile. The solution to this energy paradox is found 
in an article by C. M. De Witt and B. De Witt from 1964. The main point is that the equation of 
motion of a radiating charge is modified in curved spacetime. In the present article we explain the 
physics behind this modification, and use the generalized equation to solve the perpetuum mobile 
paradox. 



I. INTRODUCTION 

A surprising possibility of a perpetuum mobile seems 
to exist according to the usual equations of classical elec- 
tromagnetism. A situation leading to this possibility was 
described recently by Chiao He considered two ob- 
jects, one neutral and one charged, orbiting in free fall 
around the Earth, and wrote: 

The charged object will gradually spiral in to- 
wards the Earth, since it is undergoing con- 
stant centripetal acceleration in uniform cir- 
cular motion, and will thereby in principle 
lose energy due to the emission of electro- 
magnetic radiation at a rate determined by 
Larmor's radiation-power formula [using here 
SI- units] : 

q 2 

Pem =a q a 2 a q =- ^ (1) 

where Pem is the total amount of power 
emitted in electromagnetic radiation by the 
charged object with charge q undergoing cen- 
tripetal acceleration a. The energy escaping 
to infinity in the form of electromagnetic radi- 
ation emitted by the orbiting charged object 
must come from its gravitational potential en- 
ergy (which is related by the virial theorem 
to its kinetic energy), and therefore this ob- 
ject will gradually spiral inwards towards the 
surface of the Earth. 

Although everything Chiao says is not only perfectly 
reasonable, but must also be correct due to the princi- 
ple of conservation of energy, the Lorentz-Abraham-Dirac 
equation of motion (the LAD-equation) has a strange 
consequence as applied to the situation he describes. The 
covariant form of this equation is : 

f^ + V^u" (2) 

where the dot denotes covariant differentiation with re- 
spect to the proper time of the charge, and 

I* =a q (V - \a%uA (3) 



is the Abraham four-force, also called the field reaction 
force. Here a@ is the four-acceleration of the charged 
body. 

The four-acceleration vanishes for geodesic motion. 
Hence there is no field reaction force in this case. This 
means that according to the LAD-equation the neutral 
and the charged particle motions are identical, although 
the charged particle radiates, and the neutral particle 
does not. 

Since the particles return to the same point of the path 
after one period, they can continue this motion for an un- 
limited amount of time. And the charged particle would 
continue radiating all the time. Hence, this is a perfect 
perpetuum mobile. One gets radiation energy for nothing 
from a system working cyclically. 

Obviously this in conflict with the principle of conser- 
vation of energy, so something must be wrong. There 
are two equations involved in this energy paradox, the 
Larmor equation and the LAD-equation. 

One solution of this energy paradox would be that a 
charge falling freely does not radiate. One could argue for 
this by noting that a freely moving particle defines a local 
inertial reference frame, i.e. a frame in which Newton's 
1st law is valid. Let the neutral object mentioned by 
Chiao be an observer. As measured by this observer the 
charge is at rest in an inertial frame, and hence it does not 
radiate. However, a thorough analysis proves that this is 
not the solution to the paradox [H, [H, H, 0] ■ It has been 
shown that, as measured by an observer that is not falling 
freely, a freely falling charge radiates with a power given 
by Larmor's formula, and it extension to a non-inertial 
reference frame. No generalization to curved spacetime 
has been made, but since the above result shows that it 
does not matter for the radiated effect whether a charge 
is accelerated by gravitational or normal forces in flat 
spacetime, we will assume that this also holds, at least 
in the Newtonian limit of weak field and low relative ve- 
locities, in curved spacetime. 

The other equation is the equation of motion of the 
charged, radiating particle, the LAD-equation. Maybe 
there is something wrong with this. Could there be a non- 
vanishing field reaction force on a freely falling charge? 
An argument against this possibility is the following. Ac- 
cording to the principle of equivalence no local measure- 
ment should be able to distinguish between being at rest 
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in flat space and being freely falling in curved space. As 
observed by the comoving neutral object there is no force 
acting upon the charged particle. Hence it should remain 
at rest relative to the neutral object, and this is a refer- 
ence independent result. 

This argument has, however, a weak point. The prin- 
ciple of equivalence has a local character. The mentioned 
equivalence is only valid as far as the measurements does 
not reveal a possible curvature of space. So, if there is 
a non-local interaction between the charge and the cur- 
vature of spacetime, due to the non-local character of its 
electromagnetic field, this may modify the equation of 
motion of the charge. The necessity and nature of such 
a modification will be explained in the next sections. 



II. THE FIELD REACTION FORCE IN 
CURVED SPACE 



Before we can understand the solution to the energy 
paradox introduced above, we need to know the cor- 
rect form of the equation of motion of a charged par- 
ticle in curved space. It was, in fact, deduced more 
than forty years ago [1, but is not well known out- 
side a rather small group of physicists having performed 
research within this field. The reason it not that this 
knowledge is not important - the equation of motion of 
a charged particle in curved space is one of the funda- 
mental equations in classical electrodynamics - but that 
the equation is mathematically rather complicated. Nev- 
ertheless, the physics behind the equation can be under- 
stood, as we shall explain in this article. 



A. An introduction to the electromagnetic self 
force problem in curved spacetime 

Let us consider the motion of a classical charged point 
particle due to the electromagnetic force. The path of 
the particle is x'(r), and we shall also need a label for 
a general point in space-time, which will be x. To make 
the notation compact, it is usual to associate the prime 
also with components of vectors taken at a given point. 
Thus, j M is understood to be j M (x), while j M = j 1 * (x'). 

In general, the force on a body due to the electromag- 
netic interaction is given by the Lorentz equation 



dV 



=F» u f 



(4) 



where -4y is the force per unit volume on the body, the 
tensors -F M „ = A v ^ — A^^ and j 11 are the electromag- 
netic field and current density, respectively, and is 
the electromagnetic potential. The force on the whole 
body, then, is 



F»„fdV 



(5) 



In the case of a point particle with charge q four- velocity 
u and current density j M (x) = q J u^S(x, x'(r))dr, where 
r is the proper time of the particle, the Lorentz equation 
reduces to 



(6) 



The <5(x, x') = <7~'(5(x — x') used in the expression for 
the four-current is the invariant delta-function, which is 
defined such that any volume integral over x is 1 if it con- 
tains x', and otherwise, furthermore g is the absolute 
value of the determinant of the metric. 

Since the particle is charged, it not only reacts to F^ Vl 
but also acts as a source of it, according to the electro- 
magnetic wave function 



A^ a a - RitA a 



(7) 



with i? M „ being the Ricci curvature tensor. The solution 
of this equation depends on the motion of the particle, 
encoded in which is again determined by the Lorentz 
equation, so these two equations are coupled, the cou- 
pling being the particle's reaction to its own field. 



B. Green functions 



Since these equations are linear, it is possible to split 
j 1 * into points, solve the equation for each point, and then 
take the sum of the results. Furthermore, we do not need 
to know the amplitude of j M at each of these points when 
we solve the equation, as the linearity makes it possible to 
multiply by this afterwards. This way, much of the work 
in solving the equation can be done without involving 
j M itself. This is called the Green function approach to 
solving differential equations. 

In general, given some linear differential equation 



L(</,(x)) =/(x) 



(8) 



where L is a general linear differential operator (for ex- 
ample L = — J- 



dx 2 



y^j), a Green function defined by 
L(G(x,x')) =6(x-x') 



(9) 



where 6(x — x') is the Dirac delta function, solves the 
original differential equation through 



(x) = J G(x,x')f(x')dx' + S 



(10) 



where S is a surface term that can be taken to be zero 
when the integral goes over all of spacetime. 

In many cases, it is easier to find the Green function 
than the unknown function, and so it is mainly used for 
convenience, but it also serves to make the causal connec- 
tions inherent in the differential equation explicit. From 
the definition, we see that the G(x, x') is the effect on the 
spacetime point x of the point x', and thus can be seen 
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as the degree to which the two points are causally con- 
nected (as far as the differential equation is concerned). 
For example, in flat, 3+1-dimensional spacetime, electro- 
magnetic radiation moves only on the light cone, and the 
Green function for electromagnetic radiation is therefore 
zero everywhere outside of the light cone, so in this case 
G oc 5(a). The a appearing here is commonly used in 
the theory of geodesies, and is defined as a = ^s 2 , with 
s being the geodetic interval between two points (x and 
x' in this case). It is negative for a time-like interval, 
positive for a space-like interval, and zero for a light-like 
interval. 

The Green function defined above is a scalar quantity. 
However, the equation we are considering, the electro- 
magnetic wave equation, is a vector equation. A straight- 
forward generalization of equation ^ for this case is 

L^(GV(x,x')) =<5V%-x') (11) 
A»= I GV/dx' (12) 

This equation is correct [13j as it stands, but to make it an 
invariant expression, we need to replace the delta func- 
tion with its invariant version 5(x, x'), and the integral 
needs a volume element of *J~g. Rectifying this, we arrive 
at 

L^(GV(x,x')) =<5V<Kx,x') (13) 
<P=j GVf'V^x' (14) 

Applying this to the electromagnetic wave equation, we 
get 

□GV - i?^GV - - MV<5( X > x ') (I 5 ) 
A» =J G\,f'^ (16) 

where the constant factor —fj,o has been absorbed into G 
due to convention. 



C. Generalized LAD equation for curved space 

Equation (| 1 5|) differs from equation ([7]) in one impor- 
tant aspect: It is independent of j M , and thus not cou- 
pled to the motion of the particle. This is the main reason 
for introducing Green functions. It may still be difficult 
to solve, but it can be solved once and for all, and then 
applied to specific cases with equation (TTB|l. 

Let us assume that we already know the Green function 
for the wave equation; we can then use it to find the 
electromagnetic potential: 

a»{x)= J GM^')rV)V^v 

=q [ CV(x,x'(t)K'^ (17) 



Since this is an explicit expression for A^, we can insert 
it into the Lorentz equation, resulting in 

F(t) =q 2 J dr'(G pQ ,^(x'(r),x'(r')) 

- G'V;p(x'(t),xV))RVK(t) (18) 

Here, the general point x has been replaced by the po- 
sition of the particle, x', as that is the point where we 
use the Lorentz equation. The physical interpretation of 
this equation is that each point x'(r') on the world line 
of the particle gives a contribution to the force acting 
on some chosen point x'(t), the effect of x'(r') on the 
potential, and therefore on the force, being given by the 
Green function GV (x'(r), x'(r')). 

The case r = r', where emission of and reaction to 
the field happens at the same time, stands out as both 
promising due to its locality, and daunting due to the fact 
that the infinite charge density of a point charge sets up 
an infinitely strong field at its location. 

The contribution from this case was computed by De- 
Witt and Brehme Q in 1960 and Hobbs 9] in 1968 as 

- ^a q {R" u u v c 2 + R a ^u a u^u v ) (19) 

here is the Abraham four- force. The part of the in- 
tegral (fl"8|) where the "absorption" time r is earlier than 
the emission time r', can be eliminated on grounds of 
causality, and we are left with 

r(r) =1* - ±a q {R\u v c 2 + R ap u a u u») 

+ q 2 [ dr' (G pa ^ - GV ;P ) u a '( T ')uP(T) 

J —OC 

(20) 

This is a generalization of the LAD equation we used 
earlier, which stops after the first term on the right hand 
side, and is therefore a local equation. It is on this gener- 
alized equation that DeWitt and DeWitt [lOJ base their 
calculations. 

Though we have seen how the integral term, usually 
called the tail term, originates analytically, it is still not 
obvious that it can be anything but zero. We will there- 
fore take another detour before returning to the DeWitts' 
results, in order to show how this can be. 

As we remember from the derivation of the term, it rep- 
resents the force from a field that is generated at the point 
x'(r'), and later reaches the point x'(t). The behavior 
of the electromagnetic field is encoded in G, but before 
calculating, we already have some expectations about its 
form from our experience with electromagnetism: 

1. Electromagnetic fields in vacuum should move at 
the speed of light, c. 

2. It should be impossible for a particle to catch up 
with a light-like field it has emitted. 
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From these two points, we would conclude that the only 
possible way a particle can interact with itself is when 
emission and absorption of the field happens at the same 
point, so there should be no nonlocal contributions, and 
therefore no tail term! 

Our intuition turns out to be wrong on both these 
points, and both of them can give contributions to the 
tail term. 



III. THE PHYSICS BEHIND THE NON-LOCAL 
FIELD REACTION IN CURVED SPACE-TIME 

A. Overtaking light 

The latter premise for the assertion that there should 
be no tail term, namely that it should be impossible for a 
particle to catch up with a light-like field it has emitted, 
is the easiest to prove wrong, as gravitational lensing 
provides a simple counterexample. 

As an extreme case of gravitational lensing, consider 
a particle suspended at a constant Schwarzschild coor- 
dinate of r = |i?5, i.e. at the innermost stable or- 
bit. It then emits a photon perpendicularly to the di- 
rection towards the black hole. The photon will then 
reencounter the particle after one orbit. The particle has, 
effectively, caught up with the photon by taking another 
path through space-time. 

Less extreme examples of lensing also make the same 
effect possible, though a relative velocity between the 
mass and the particle will be necessary in those cases. 



B. Slow light 

In the above case, it was possible to catch up to 
light because light took another path through spacetime, 
though it still moved at the speed of light. It is more 
difficult to see how light can end up moving more slowly 
than c through empty space. However, this turns out to 
be the most important contribution to the tail term. 

To see how the wave equation can lead to parts of the 
field moving more slowly than the speed of light, and 
thus causing dispersion, so that an initially sharp pulse 
develops a trailing tail, let us first consider the spherically 
symmetric case in flat spacetime, and vary the number 
of spatial dimensions. These cases provide one of the 
conceptually simplest examples of dispersion, and will 
turn out to be mathematically related to cases in curved, 
3+1-dimcnsional spacetime. 

In flat spacetime, it is possible to decompose the wave 
equation, which is a vector equation, into one scalar wave 
equation for each dimension of space-time, with all the 
scalar equations being independent. These scalar equa- 
tions have the form 

<P„ - -P (21) 



and solving this equation is equivalent to solving the vec- 
tor equation. In curved spacetime, this is no longer the 
case, but the solutions to the vector and scalar equations 
generally display similar behavior. 

In polar coordinates with spherical symmetry, equation 

© takes the form (-^^ + £t + 4> = ~Pr, 

with p r being the angular integral of p, and d being 
the number of spatial dimensions. The equation for the 
Green function corresponding to this is 

/ Id 2 d 2 d - 1 d \ 

\ c 2 dt 2 dr 2 ^~ r dr ) ()()() 

Note that the only thing separating the cases with dif- 
ferent number of dimensions, is the factor in front of the 
term, but this factor is critical for its behavior. It 

or ' 

possesses simple solutions for the cases d = 1 and d = 3, 
but is complicated to solve for other cases, especially the 
even ones. Brown fill ] discusses this thoroughly in an 
article on the Huygens' Principle. He only considers the 
homogeneous version of the equation, 

fid 2 d 2 d-1 d\ , „ 

but since the right hand side of equation ([22]) is given by 
delta functions, they only differ at the initial moment, 
and so the solution to equation (|2"2"|) can be found by 
solving equation f|23|) with an appropriate initial condi- 
tion. 

To see that the case for 1 spatial dimension and our 
familiar 3 spatial dimensions are exceptionally simple, 
Brown suggests that one makes a change of variable by 

d— 1 

introducing ip — r ? ^ (our notation switches the role of 
4> and ip compared to his), which leads to the differential 
equation 

(~^dt 2 ~lfr 2+ 4^2 a^)^" ^ 

The first order term vanishes for d = 1 and d = 3, giving 

1 ^ 2 / ft 2 1 

the simple equation -g^f- = , which can be rewritten 

as (f - hit) (£ + \i) V> = °> and 80 obviously has 
solutions 

ip = f(r - ct) + g(r + ct) (25) 

for any function / and g [l4| and then It is easy to see that 
this solution consists of signals moving at a single, well- 
defined speed: c. This exceptional case (d = 3) is what 
our intuition for the wave equation, and thus behavior 
of light, is based on, and the behavior of electromagnetic 
fields in other numbers of dimensions therefore end up 
being counterintuitive. For a general d, we cannot expect 
the first order term to cancel, and therefore do not end 
up with an an equation that leads to this kind of simple 
propagating solution. 

Brown goes on to consider the case d = 2, and by as- 
suming separability, <fi(r,t) = f(r)g(t), rewrites equation 
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Figure 1: Time series for the solution of the equation 
{-^dt +dr + ^-d r )G = -2cS(t)S(r) with d = 2 and 
Gaussians instead of delta functions, and c — 1. r runs along 
the horizontal axis, and G along the vertical one. This is a set 
of superimposed curves, one for each time step. The general 
behavior is that of a pulse propagating to the right, but the 
value behind the pulse does not reach 0, making the final re- 
sult a lopsided Gaussian. Since we have spherical symmetry, 
the time series in this figure corresponds to an initial central 
peak evolving into a spreading circular "hill" , with a nonzero 
"elevation" inside, which is the "tail". 



(l23f as 

9 d 2 f df 9 „ „ d 2 q , , . 

r ^ +r £ +rf=0 J=- ca (26) 

The equation for / is the zeroth order Bessel equation, 
and is solved by the corresponding zeroth order Bessel 
functions Jo(r) and Yq(t). By using the integral repre- 
sentation for Jo, Jo( r ) — f / °° sin (cosh(#)r) d8, and a 
solution g(t) = sin(ct) of the equation for g, Brown finds 
a solution for <fi to be 

1 f°° 

4> =— I [cos (cosh(6>)r — ct) — cos (cosh(6')r + ct)] d9 

T Jo 

(27) 

This, he points out, is a superposition of waves traveling 
with speeds of v = ^§7gjj which takes all values from 
to c. So in this case, light moves at the speed of light, as 
well as all lower velocities*. 

To see how this applies to the actual Green equation, 
which differs from the equation considered by Brown by 
being inhomogeneous, we present numerical solutions for 
the cases of 2, 3 and 4 spatial dimensions, and with delta 
functions approximated by Gaussians. The result can be 
seen in the figures [TJH 

From these solutions, we see that the factor in front 
of the -ftp term of equation (|2"2")l controls how fast the 
Green function falls after the wave front has passed. For 
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Figure 2: As figure [T] but for d = 3. This time the value be- 
hind the pulse does reach 0, and we get a normal propagating 
Gaussian, corresponding to a propagating delta delta func- 
tion. Since the amplitude is zero before and after the wave 
front, radiation propagates only on the light cone in this case. 
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Figure 3: As figure [TJ but for d = 4. This bears more resem- 
blance to the case of n = 2 than n = 3, in that the amplitude 
behind the wave front does not reach zero, but this time it is 
due to it falling too rapidly, overshooting zero. 



d = 2, which gives a factor of the fall after the wave 
front is less than the rise that preceded it, so the function 
develops a tail. For d = 4, with a factor of -, the fall is 
too great, and the solution overshoots zero, again leaving 
a tail. The case d = 3 has the exact balancing between 
fall and rise needed to avoid a tail. We again see that not 
having a tail is not a general feature of the wave equation, 
but a special case. 

To see how this applies to curved 3+1-dimensional 
spacetime, we now consider the form of the spherically 
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symmetric scalar wave equation in the Schwarzschild 
spacetimc: 



£J?1 

c 2 dt 2 



1&_ 

(3 dr 2 



1+0 d 
r dr 



G = - 2cS(t)S(r) (28) 



with the shorthand notation (3 = 1 — —. This equation 
is expressed in terms of Schwarzschild radial coordinates, 
which behave as flat-space polar coordinates only in the 
limit r — > oo. This prevents us from directly comparing 
it to equation (122p . To remedy this, one can introduce 

orthonormal coordinates r = ^/l — —t, p = -, r 

takes the form 



around r = r . With these, equation 



i d 2 

c 2 dr 2 



dp 2 



2-3 



Rs 



R% 

Tf d 



dp 



G 



2cS{t)S( p ) 
(29) 



This equation is of the same form as equation (|22j) . but 
this time the factor in front of the first order radial deriva- 
tive varies between (in the limit — > Rs) and 2 (in the 
limit ro — > oo). Except for in the latter limit, these fac- 
tors are too small to make the Green function reach zero 
behind the wave front. Hence, there will be dispersion. 

The right hand side of this equation is problematic, as 
the source term is in the middle of the black hole. This 
does not matter much, though, since we are interested in 
the fate of an initially sharp wave front, and we only need 
the homogeneous version of the equation to find the time 
evolution of such a wave. One could, for example, start 
with an initial wave concentrated in a thin shell centered 
on the black hole (a shell because we assume spherical 
symmetry here). The homogeneous equation then tells 

us that this will evolve as in a 3(1 — — ) + ^-dimensional 
spacetime, and thus disperse accordingly. 

Though we considered the scalar wave equation here, 
one will in general find a tail for the Green function in 
curved spacetime also for the electromagnetic wave equa- 
tion 0. 

The existence of a tail of the Green function means that 
parts of the field spread at speeds slower than the speed 
of light, and thus move inside the light cone. In general, 
parts of the field will be found to move at all speeds slower 
than the speed of light [H, pg. 1 1] . When a sharp wave 
front enters an area of curvature from flat spacetime, it 
will begin developing a trailing tail stretching not only to 
the point where the curved area starts, but also beyond 
this, since dispersion also leads to part of the field moving 
in the opposite direction to that of the initial wavefront. 
This part can propagate back to the particle and interact 
with it. 

There are, then, 4 ways for the particle to interact with 
its own field: 

• Locally, where the particle reacts to the field it just 
emitted. 



• By taking a different path through spacetime than 
the field, which lets it encounter even the parts of 
the field that move at the speed of light. 

• By catching up to parts of the field it emitted ear- 
lier, because the field is moving slower than c. 

• By having parts of the field scattered back at it. 



IV. THE SOLUTION OF THE ENERGY 
PARADOX 



The question is: How does energy conservation come 
about in a situation where a charged particle moves in 
the Schwarzschild spacetime, say in the equatorial plane, 
under the action of gravity alone? Assuming the particle 
moves in a circular orbit [la] far away at a distance r from 
a mass M, it will have a motion given by 



r = ref 
a 



v = roe; 



ru> 2 ef 



with a = |a| = GMr 2 , giving a radiated effect of 



P 



q 2 G M 2 
67rr 4 



(30) 



(31) 



by the use of equation {T]), and the particle needs to lose 
energy at the same rate for energy to be conserved. 

The answer to the question above is contained in 
the generalized LAD-equation for curved space, equation 
(|2"0"|) . In the present case the local term containing the 
Ricci curvature tensor does not contribute since this ten- 
sor vanishes in the external Schwarzschild spacetime. It 
is the non-local so-called tail-term that matters here. 

Before we discuss the solution of the paradox, let us 
note an interesting conceptual difference between the spe- 
cial relativistic version of the equation of motion as de- 
duced by Lorentz and Abraham on the one hand, and 
the fully covariant version of it deduced by Dirac. In 
the present case we are not concerned with velocity- 
dependent effects. Hence we may start by looking at 
the non-relativistic version of the equation, 



(32) 



The special relativistic (and Newtonian) acceleration is 
absolute and arises when the particle is acted upon by 
an external force (or during pre-acceleration or run-away 
motion) . Hence, during circular motion around the Earth 
under the action of the gravitational force, there is an ex- 
ternal gravitational force F cxt = G1 ^i m acting in the ra- 
dial direction. The change of the acceleration is directed 
opposite to the motion, and so the field reaction force 
acts like a frictional force, decreasing the velocity of the 
charge, which thereby starts spiraling inwards. The radi- 
ated energy is accounted for by a decrease in mechanical 
energy. 
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As explained by DeWitt and DeWitt [T(| , the contents 
of the covariant equation ^ interpreted within the gen- 
eral theory of relativity is different. According to this 
theory, gravity is not reckoned as a force. In Newton's 
theory and in special relativity we say that a freely falling 
particle is acted upon by a gravitational force. In Ein- 
stein's general theory of relativity we say that a freely 
falling particle is not acted upon by any force. Such a 
particle has vanishing four-acceleration. DeWitt and De- 
Witt then say: 

When a charged particle is accelerated by 
means of nongravitational forces, the electric 
field lines which emanate from it bend and 
redistribute themselves in the vicinity of the 
particle (i.e., within a distance of the order of 
the classical radius) in such a way as to ex- 
ert, on the average, a net retarding force, over 
and above the force of inertial reaction. With 
purely gravitational forces, however, this is 
not the way things happen. The field in 
the immediate vicinity of the particle tends 
to fall freely with the particle, and although 
it suffers a local tidal distortion characteris- 
tic of an explicit occurrence of the Riemann 
tensor . . . the net retarding force due to this 
distortion is zero when integrated over solid 
angle [l6j. 

Hence, the usual Abraham force plus local contributions 
to the field reaction force from terms depending upon the 
curvature vanish for free particles. However, DeWitt and 
DeWitt continue, saying: 

The deviation of the particle motion from 
geodesic when F^ t — is caused not by the 
local field of the particle but by a field which 
originates well outside the classical radius and 
which is manifested by the nonlocal term of 
equation [ ([20| ], 

They calculated this by considering the whole electro- 
magnetic field (and not just the part of it in the neigh- 
borhood of the particle, as one normally would) using 
the weak field approximation and the assumption of a 
static gravitational field. The result is that the curved 
spacetime around a massive object disperses and scatters 
the electromagnetic field. Thus, it is possible for the field 
to emanate from the particle, propagate until it reaches 
an area with non-negligible curvature, be scattered there, 
and then propagate back to the particle. The DeWitts 
write 

Physically [this] arises from a back-scatter 
process in which the Coulomb field of the par- 
ticle, as it sweeps over the "bumps" in space- 
time, receives "jolts" which are propagated 
back to the particle. 

This non-local interaction between the particle and its 
field provides a mechanism for reducing the particle's en- 
ergy and balancing the energy budget. But as opposed 



to the radiation reaction for nongravitational forces, this 
force is determined not by the particle's current motion, 
but by its motion at some earlier time, and so there is no 
reason to expect these to match, or even resemble each 
other. 

This non-locality also makes the equation of motion a 
time delay differential equation, a difficult class of dif- 
ferential equations that often lack analytical solutions. 
However, by going to the non-relativistic limit, DeWitt 
and DeWitt succeed in evaluating the extended field re- 
action for a point charge in free fall near some point mass, 
and near some arbitrary mass distribution, the result be- 
ing 

F =Fc + Fjvc 
F c =q 2 GM^ 

2 / 1 rr\ 

F NC = - r 2 GM ^ - 3^J • r (33) 

where M is the mass of the source of the gravitational 
field and q is the charge and r is the separation 3- vector 
between the charge and the mass as measured in the rest 
frame of the mass. This can be written 

Fc = — mVtp 

Fnc = ~ \q 2 i ■ VV0 = -^ 2 ^-(V0) 
3 6 at 

2 m J (r — r'y 

m =- G J |7=^ V (34) 

where p is the mass density and m<p is the normal grav- 
itational field. The mathematical expressions for Fjvc 
are explained in Appendix [A] DeWitt and DeWitt write 

Fnc is the nonconservative force which gives 
rise to radiation damping. Owing to its de- 
pendence on velocity it is small in magnitude 
compared to the force Fc which is conserva- 
tive. ... It [the conservative part] arises from 
the fact that the total mass of the particle is 
not concentrated at a point but is partly dis- 
tributed as electric field energy in the space 
around the particle. 

So F^vc is the force that should correspond to the ra- 
diated effect in equation ([31) . However, this effect was 
calculated under the assumption of a circular orbit, which 
is not possible under the influence of a nonconservative 
damping force. But if Fnc is very small, as it should be, 
because the field has traversed a distance of 2r before 
returning to the particle, and only a small part of the 
field is scattered back in the first place, then the parti- 
cle motion will be dominated by the force of gravity, so 
■Ftotai ~ Fq = -mV4>. 

By inserting this into the expression |3"4"1) for F^c, we 
see that the nonconservative force, which with its de- 
pendence on velocity and the local shape of <fi is quite 
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different from the Lorentz- Abraham force Fla = a q k, 
still reduces to it in the non-relativistic limit. The work 
this force does on the particle is given by P = F^c 1 v > 
and using equation ([3D)) , we find 



P =a q k • 



2 4 

-a„r u) 



q 2 G 2 M 2 
67rr 4 eoC 3 



(35) 



which is of the same size but opposite magnitude as the 
radiated effect in equation (|31[) , meaning that the particle 
loses energy at the same rate as it is radiated away, and 
energy is thus balanced. 

This applies not only for the circular orbit considered 
here; DeWitt and DeWitt showed (still using the as- 
sumption that the motion is dominated by the effects 
of gravity) that the average energy loss rate during any 
closed orbit in free fall matches that given by the Lar- 
mor formula. This follows from the fact that the non- 
local, nonconservative, backscattering-based force F N q 
reduces to normal Lorentz- Abraham force Fla in the 
non-relativistic limit. This, then, is what balances the 
energy budget, and prevents an electron in orbit around 
the Earth from being a perpetual source of energy. 
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Appendix A 



Given two vectors and two obvious products 
can be formed: the inner product P = A^B^, and the 
outer product Q^ v = A^B V . DeWitt and DeWitt indi- 
cate the former with the • operation, and the latter by 



adjoining two vectors without any intervening notation. 
Thus, the rr seen in equation (j3"3"|) is 



n 3 



T\r\ r x r 2 rir 3 
r 2 r\ r 2 r 2 r 2 r 3 
,r 3 ri r 3 r 2 r 3 r 3 . 



(Al) 



and (rr)-r = r(r-r). The matrices in equation (f3"3"| simply 
serve to let one write the r as a single factor outside the 
parentheses; otherwise, it could be written without them 
as 



-q GM ( -ir - 3r— jp 



(A2) 



The above also applies to the vector operator V, which 
is (d x ,d y ,d z ) in Cartesian coordinates. The inner prod- 
uct formed from V is V 2 = V • V = d 2 + d 2 + d 2 , but 
one can also form the outer product 



r d x d x d x d y d x d z 

VV = ( 8yd X dydy 3 y3 Z 

,d z d x d z d„ d z d z 



(A3) 



While V 2 = in areas where the mass density p = 0, 
the same does not apply to VV0; otherwise, F^c would 
be at the particle's position, and we would have no 
extended field reaction. 

Note also that the fact that the gravitational field is 
static has been used in the transition v • VV</> = JjV0. 
In general, the total time derivative of some field £ as 
experienced by a particle moving with velocity v is given 
by 



— P=— f 

dr~dr~ 

but the partial derivative of £ 
is here. 



V£ 



(A4) 



with respect to time 
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